Deep-Elastic pp Scattering at LHC from Low-x Gluons 

M. M. Islam"-, J. Kaspar^^'^ and R. J. Luddy'^ 

"Department of Physics, University of Connecticut, Storrs, CT 06269, USA, (islam@phys.uconn.edu) 
*" TOTEM Collaboration, CERN, Geneva, Switzerland (Jan.Kaspar@cern.ch) 
"^ Institute of Physics, Academy of Sciences of the Czech Republic, Prague 
Department of Physics, University of Connecticut, Storrs, CT 06269, USA, (RJLuddy@phys.uconn.edu) 

Abstract 

^^ Deep-elastic pp scattering at cm. energy 14 TeV at LHC in the momentum transfer range 4 GeV^ ^ |t| ^ 
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10 GeV is planned to be measured by the TOTEM group. We study this process in a model where the 



^vg , deep-elastic scattering is due to a single hard collision of a valence quark from one proton with a valence 

quark from the other proton. The hard collision originates from the low-x gluon cloud around one 
f*> I valence quark interacting with that of the other. The low-x gluon cloud can be identified as color glass 

■^ ' condensate and has size ~ 0.3 F. Our prediction is that pp da/dt in the large |i| region decreases 

smoothly as momentum transfer increases. This is in contrast to the prediction of pp da/dt with visible 
oscillations and smaller cross sections by a large number of other models. 
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p i' A major discovery from deep-inelastic ep scattering at HERA is the dramatic increase of 

^ . low-x gluon density in a proton [1]. The HERA data have been well described by the saturation 

model of Golec-Biernat and Wiisthoff (GBW model) [2,3]- The prediction of geometric scaling 

in this model for low-x region (x < 0.01) has been impressively confirmed by Stasto et al.[4]. In 

^ . GBW model, the virtual photon splits into a dipole (a qq pair) and the low-x gluons of a proton 

lO ! present an effective black disk to the dipole in the transverse 6-plane. In this paper, we present 

}£? ' results of our investigation of deep-elastic pp scattering at LHC using key premises of GBW 

Q ■ model pertaining to the low-x gluons. 

Tj- ■ pp elastic scattering at LHC at center-of-mass energy 14 TeV and momentum transfer 

O , range |t| ~ - 10 GeV^ is planned to be measured by the TOTEM group [5,6]. By deep-elastic 

Q I pp scattering — we mean elastic da/dt distribution in the momentum transfer range 4 GeV^ < 

|t| < 10 GeV^. In this range, a proton probes the other proton at a transverse distance b ~ q^^ 
(Q = VW) ^11^ therefore at distances of one tenth of a fermi (F) or smaller. High energy pp 
elastic scattering has been studied by us in a model of proton where the proton has an outer cloud 
C^ ■ of qq condensed ground state (size ~ 0.86 F), an iiuier shell of baryonic charge density (size 

~ 0.44 F) and an inner core of valence quarks (size ~ 0.2 F) (Fig. 1) [7]. Large \t\ pp elastic 
scattering in this model is due to a hard collision between a valence quark from one proton with 
a valence quark from the other proton — where the collision carries off the whole momentum 
transfer \t\ (Fig. 2). In a previous paper, we have studied the case where this collision is due to a 
hard pomeron (a BFKL pomeron plus its next-to-leading order corrections) [8]. We now want to 
study the case where the low-x gluon cloud of one valence quark interacts with that of the other 
valence quark leading to the deep-elastic scattering. 

The valence quarks in our model are color neutral quasi-quarks, because of the gluon 
clouds around them. To describe the valence quark-quark (q-q) scattering, we consider first an 
eikonal description for a purely absorptive process where the scattering amplitude is given by 
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f{s,t)=ipW^JdH e''-^[l -exp(-xi(s,6))]; (1) 

Xi(s, b) is the imaginary part of the eikonal and the hat C) refers to valence quark-quark scatter- 
ing. From Eq.(l), we obtain the Born (single scattering) amplitude 

fB{s,t) = ipW ^ J dH e^'-^ Xi{s,b).^ (2) 

Applying the optical theorem: Im T{s, 0) = ^^ (Ttotls) to Eq. (1) and, in Born approximation, 
to Eq. (2), we find 

atot(s) = 2jdH[l- exp(-xi(s, b))] (3) 

and 

aB{s) = 2jdHxi{s,b), (4) 

where cjb(s) is the total cross section in Born approximation. Let us next consider a factorizable 
eikonal: xi{^,b) = g{s) ^ F{b), where g{s) is an unknown real function of s and F{b) is a 
normalized 6-profile: / d^6 F{b) = 1. Then Eq. (4) leads to 

^bIs) = g{s), (5) 

i.e. the unknown function g{s) has to be the total cross section in Born approximation. The 
corresponding total and Born amplitudes from Eqs. (1) and (2) are 

fis,t)=ipW^JdH e'b-g^l - exp{-Ms)i Fib))] (6) 

and 

fBis,t) = ipW&Bis) i^ J dH e^'-^Fib). (7) 

Eqs. (6) and (7) then show that to obtain elastic valence quark-quark scattering amplitudes (multi- 
ple or single), we need two physical quantities: asis) and the 6-profile of valence q-q scattering. 

We view each valence quark as having a low-x gluon cloud around it described by a gluon 
density g{x, Q'^{x)) and a transverse space profile f{b). (Our choice of the scale Q'^{x) in gluon 
density comes from the GBW model, as will be seen later.) By noting that (?(x, (5^(x))da; f{b)d^b 
is the number of low-x gluons in the gluon cloud of a valence quark lying between fractional mo- 
mentum X and X + dx and in the 6-plane area d^b, and by introducing a low-x gluon-gluon 
production cross section, we can write down a simple expression for valence q-q total cross 
section (which we take as Born cross section): 

<7b(S) = jdH j^^ dxi g{xuQKxi)) JdHi f{bi)J^^ dx2 g{x2, QK^2)) 

xIdH2f{b2)d{xi X2s)6(^\b + bi -h); (8) 

here (j{xi X2 s) is the low-x gluon-gluon production cross section at cm. energy squared xiX2S 
(s is the cm. energy of the colliding valence quarks: s = {p + k)"^ ~ p+k-; see Fig. 2). We 
refer to gluons lying in the range < x < Xc (Xc = 0.01) as low-x gluons. (7b(s) can now be 
written as 

aB{s) = agg{s) JdHF{b), (9) 

where 

'^gg(«") = lo' dxj g{xi , Q^{xi)) Jo" dx2 g{x2, QK^s)) ^{xi xg s) (10) 

and 

F{b) = j dHi f{bi) j dH2 f{b2) 5^^\b +bi -b2). (11) 

Equation (10) provides a low-x gluon-gluon total cross section and Eq. (11) provides the corre- 
sponding 6-profile F{b). Comparing Eq. (9) with Eq. (4), we notice that 



2xi(s,6) = agg(s)F(6). (12) 

Furthermore, the low-x gluon 6-profile f{b) of each valence quark represents probability den- 
sity in the 6-plane and therefore is normalized: J dj^b f{b) = 1. Equation (11) then leads to 
/d^6 F{b) = 1 and from Eq. (9), we obtain cjb(s) = o'gg(s)- Finally, from Eqs. (6) and (7), 
we find the q-q scattering amplitude to be 

f{s,t)=ipW^!dH e^Hl -exp(-|agg(s) F{b))] (13) 

and the Born amplitude to be 

fBis,t)=ipWa^^{s)^JdH e^'-'lFib). (14) 

At this point, let us recall that our goal is to study deep-elastic pp scattering due to a single 
hard collision of a valence quark from one proton with that of the other (Fig. 2), where the 
collision carries off all the momentum transfer and originates from the low-x gluon cloud of one 
valence quark interacting strongly with that of the other. We, therefore, focus on evaluating the 
Born amplitude, Eq. (14). The expression for o-gg(s), Eq. (10), shows that to evaluate it, we need 
to know the x-dependence of the gluon density g{x, Q'^{x)). The expression for F{b), Eq. (11), 
shows that for its evaluation we need to know the 6-profile f{b) of the low-x gluon cloud of a 
valence quark. 

We now want to point out that the analyses of deep-inelastic ep scattering in the GBW 
model, in fact, provide precise information on the x-dependence of the gluon density g{x, Q'^{x)) 
and also insight into the 6-dependence of the profile f{b). To understand this connection of GBW 
model with our work, we note that in deep-inelastic ep scattering, the virtual photon with large 
Q^ probes very small distances inside the proton, and therefore the low-x gluon density it finds 
is that due to the valence quarks - given our model of the proton structure where valence quarks 
form the core of the proton (Fig. 1). Furthermore, since the valence quarks are quasi-quarks, and 
as such are quite independent of each other, the low-x gluon density gp{x, Qg{x)) of the proton 
probed by the photon can be considered as the sum of the low-x gluon densities of the valence 
quarks, i.e. 

g,{x,Q^,{x))=3g{x,Q^,{x)). (15) 

In the GBW model, the virtual photon converts into a small qq dipole which traverses the 
low-x gluon cloud of the proton, and the dipole-proton total cross section in impact parameter 
representation is 

atr'"ix,r^) = 2jdH[l- exp(-r2 ^as x g,ix,fj.^) | f{b))], (16) 

where f{b) represents the gluon distribution in the 6-plane (J d'^b f{b) = 1). GBW next intro- 
duce a saturation scale Q'^{x) related with the gluon distribution: 

4^a,x5p(x,Q2(a;)) = Q2(^) (17) 

and take 

QKx) = Q^,i^)\ (18) 

where Xc is the cut-off value of x mentioned earlier. GBW take A = 0.29 from their analyses 
of data. Equations (17) and (18) show that as x ^ 0, x gp{x, Ql{x)) increases rapidly as x~^. 
{Qlx^ in Eq. (18) is connected with Ql x^ of GBW: Ql x^ = Ql x^). From Eqs. (15), (17) 
and (18), we find 



xg{x,Ql{x)) = C{^)\ (19) 

where C is a constant. Inserting this in Eq. (10), we obtain 

^99{s) = C^ Jt ^(fl)' St ^mf^{x, X, s). (20) 

This equation shows that the dominant contribution on the right-hand side comes from the region 
xi ^ 0, a;2 ^ 0; that is, from the threshhold region of (t(xi X2 s). Equation (20) can be written 
as a single integral over the cm. energy squared C, of the gluon-gluon collision: 
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99i^) = c' lo'' fc-frm- (21) 

As the dominant contribution to the integral comes from the threshhold region, we can approxi- 
mate a{() ~ o"(Co) ^(C ~ Co)» where (q is the threshold value and fj(Co) is the threshold gluon- 
gluon cross section. The above equation then becomes 

a,,{s) = C' a{Co)lf f (f )Mn(f ). (22) 

The integral on the right-hand side of Eq. (22) can be carried out and leads to the result 



(Tgg{s) ~ s^ Ins for large s. We note that, if the threshold gluon-gluon cross section is dominated 
by production of dijets, then a'(Co) = crdijct(4m|), where rrij is the average mass of a single jet. 

We next turn to the 6-dependence of the GBW model. Using Eq. (17) in Eq. (16), we 
obtain 

otr'^xy) = 2jdH [1-eM-r'f QHx)ifim- (23) 

In GBW model, the low-x gluon cloud of a proton is taken as a black disk: f{b) = (vri?^)^^ 6{R— 
b), where R is the black disk radius. Identifying do as the black disk cross section, ctq = 2tiR^, 
Eq. (23) becomes 

a'^r'^x^) = 2jdH[l- exp(-l r2 Q2(x) 0(R - b))] 

= ao[l-exp(-lr2Q2(x))], (24) 

which is the total dipole cross section formula in GBW model. This result, however, is prob- 
lematic. When the dipole-proton cm. energy squared goes to infinity (i.e., ^ -^ co), Eq. (24) 
leads to atotix^r"^) -^ ctq (a constant). This, however, does not meet the expected saturation 
of the total cross section: (t^^° "^(x, r^) ~ In^^ as x ^ 0, based on the Froissart-Martin (FM) 
bound [9-11]. On the other hand, if f{b) ~ ^-mb ^^j. j^gg ^^ (j,gQ (-j^g p]y[ bound is saturated. 
Furthermore, numerical investigation of the QCD nonlinear evolution equation by Gotsman et al. 
shows that f{b) falls off as e^™^ for large b, when they exclude the kinematical region of large 
dipoles [12]. To address this problem, we require the low-x gluons of a valence quark to have 
a 6-plane distribution f{b), such that: 1) /(6) falls off as e^™^ for large b, 2) has finite value at 
6 = 0, and 3) satisfies the normalization condition J d'^b f{b) = 1. 

Defining the Fourier transform of f{b) by f{q): 

f{q)=JdHe^^^-^f{b), (25) 

we obtain from Eq. (14) 

fB{s,t) = ipWa^^is)^ p{q). (26) 

As we remarked earlier after Eq. (22), CTgg{s) ~ s^ Ins. However, since the power A is not 
precisely known, we drop the Ins term and write Eq. (26) as 



TB{s,t)=is^^gS^f{q), ^ (27) 

where 7gg is a real positive constant. Requiring TB(s,t) to be crossing even, i.e. high energy 
elastic valence q-q and q-q scattering to be the same, modifies Eq. (27) to 

fB{s,t)=is^^^{se-'^)^P{q). (28) 

We next choose a general form for the low-x gluon distribution in the 6-plane: 

fib) = ^r(^i^rK,imb), (^>0) (29) 

which satisfies the three conditions we require of f{b). Equation (28) then becomes 

Tb (s, t) = z s 7gg (s e"* i )^ \2(M+i) for ^w-x gluons. (30) 



As we mentioned earlier, previously we investigated valence quark-quark scattering where 
the valence quarks interacted via the hard pomeron, i.e. BFKL pomeron (reggeized gluon lad- 
ders) plus its next to leading order corrections [8]. Because of significant leading order correc- 
tions, the power of s: s'^ with lo = wbfkl is substantially reduced to u; ~ 0.13 - 0.18 as argued 
by Brodsky et al.[13]. The Born amplitude we used was 

Tb(s, t) ~ i s 7hp (s e^* 2 )'^ -? — - — ^ for the hard pomeron (31) 



{mil = "Tg , 7hp = 7qq/"^h using earlier notation). Comparing Eqs. (30) and (31), we notice 
that the structure of the two Born amplitudes are the same. This enables us to carry out present 
quantitative calculations the same way as before [8]. pp elastic scattering amplitude originating 
from the valence q-q scattering amplitude (Eq. (30)) due to low-x gluons is 

rqq(s, t) = zs 7gg (s e~^ 1)^ , ^'l\"2S+i) > (32) 



where s = {P + K)"^ is the square of the cm. energy of pp scattering (Fig. 2) and J^'^{q±, A) 
are structure factors that arise when the momentum distributions of the colliding valence quarks 
inside protons are folded in. Physically, each structure factor stems from the confinement of 
valence quarks inside the proton. 

Our predicted pp da/dt at LHC at ^/s = 14 TeV for the whole momentum transfer range 
|t| = - 10 GeV^ due to the combined diffraction, a;-exchange, and single hard q-q scattering 
from low-x gluons (i.e. q-q scattering due to the Born amplitude: Eq. (30)) is shown by the solid 
line in Fig. 3. Also shown in the figure are da/dt due to each of these three processes separately. 
As before, the parameters of the model were required to provide good fits of pf3 da/dt at ^/s 
= 546 GeV, 630 GeV and 1.8 TeV as well as satisfactory descriptions of cJtot(s) and p{s) as a 
function of ^/s. In our present calculations, the parameters m and fi of Eq. (21) have the values 
m = 1.67 GeV, and fi = ■^, while the parameter A responsible for the sharp rise of low-x gluon 
density has the value A = 0.29 as given by GBW. 

High energy pp scattering originating from valence quark-quark scattering has also been 
discussed by Kovner and Wiedemann [14] in a model where the proton is taken as a loosely bound 
state of three valence quarks. At some initial high energy, they view each valence quark as a 
black (i.e. totally absorbing) disk with a gray periphery. The black disk originates from saturated 
gluon density in the central region of a valence quark, while the gray periphery is due to the long 



range Coulomb-like field of massless gluons. In valence q-q scattering, the interaction between 
the two gray regions leads to a power like growth of the total cross section and in their model 
explains the pre-asymptotic behavior of atot — 21.70 s^'^S'^s due to the soft pomeron [15]. With 
increasing s, the black central region of each valence quark grows in transverse size. Comparing 
the valence q-q scattering in KW model with that of ours, we notice that in our model it is the 
interaction of the low-x gluons of a valence quark with those of the other in the 6-plane (Eq. (8)), 
which leads to a power like growth of the Bom cross section <Jgg{s) ~ s'^; here A measures the 
rapid growth of the gluon density: x g{x,Ql{x)) ~ x~^. When unitarized using the eikonal 
description (Eq. (13)), the power growth of the cross section will be tamed at asymptotic s by 
the nonperturbative exponential fall-off of F{h) ~ e^™ ^ and will lead to a black disk of radius 
^BD — (A/m) In (s/so) and crtot(s) — 2 tt [{X/m) In (s/sq)]^. Furthermore, the black disk 
will have an outer gray boundary region where the profile function T{s, h) = \ — exp(— xi(s, h)) 
falls from r(s,6) = 1 (6 < 6bd) to r(s,6) = (6 > 6bd)- As s increases, the black disk 
radius will grow in transverse size and the gray region will move outward. Since f{h) d^6 is the 
probability of a gluon to be in an area d^ 6 at a distance b from the center of the valence quark 
(Eq. (8)), we can calculate the average size of the low-x gluon cloud surrounding a valence quark 
in the transverse plane; namely, < 6^ >^/^ where < 6^ >= /6^ f{h) d^h. Using the values of 
the parameters m = 1.67 GeV and ^ = ^ for f{h) (Eq. (20)), we obtain < I? >i/2= o.27 F, 
somewhat larger than the 0.2 F size used by KW in their estimate of atot, but in agreement with 
the value 0.3 F suggested by Kopeliovich et al. [16]. 

High energy pp, pp elastic scattering have been quantitatively studied by Block et al. in 
a QCD-inspired phenomenological model [17]. In the cm. system, a proton is viewed in their 
model as a disk in the 6-plane with a gluon distribution and a valence quark distribution. Both 
6-plane distributions are described by dipole form factors. Block et al. envisage gluon-gluon, 
gluon-quark and quark-quark interactions in elastic scattering, and derive the full scattering am- 
plitude using an eikonal description. Their work has been followed up by Luna et al. [18] who 
introduce a dynamical gluon mass as an infrared mass scale. Calculation of the Born cross sec- 
tion in these models is similar to ours (Eqs. (9) and (20)), since these models incorporate the 
same low-x behavior of the gluon density. However, these models do not associate low-x gluons 
with the gluon cloud of a valence quark. Instead, the low-x gluons are taken to belong to the 
whole proton. 

We note that the low-x gluons we are considering are specified at the saturation scale 
Ql{x) (Eq. (17)) where Ql{x) = Q^{xc/x)'^ (Eq. (18)). In the phase space plot: ln(l/x) versus 
InQ^, gluons lying in the region bounded by the saturation line InQ^ = lnQ^{x) and the line 
InQ^ = lnQ^{x) (separating the nonperturbative regime from the perturbative regime) are re- 
ferred to as forming Color Glass Condensate (CGC) [19]. The low-x gluons in our calculations, 
therefore, lie on the boundary of the CGC region. The CGC region, however, can be extended fur- 
ther in InQ^ using geometric scaling; namely, Q'^ < Q^ci^), Q^ci^) = Qs{x)[Qs{x)/ Q^]^-^'^ 
[19]. As a consequence, the low-x gluons we are considering actually lie well within the CGC 
region and the low-x gluon cloud around a valence quark can be identified as CGC. Deep-elastic 
valence quark-quark scattering is then due to CGC around one valence quark strongly interacting 
with that of the other. 

Our predicted pp elastic da/dt at LHC at cm. energy 14 TeV and |t| = - 10 GeV^ 



shown in Fig. 3 reflects the structure of the proton given in Fig. 1. In the small |t| region, the 
outer cloud of one proton interacts with that of the other giving rise to diffraction scattering. In 
the intermediate |t| region (1 GeV^ ^ |t| ^ 4 GeV^), one proton probes the baryonic charge 
density of the other proton via vector meson lo exchange. In the large |t| region (4 GeV^ ^ |t| ^ 
10 GeV^), one proton scatters off the other proton via a single hard valence quark-quark collision 
(Fig. 2). The valence quark-quark collision underlying Fig. 3 originates from the low-x gluon 
cloud of one valence quark strongly interacting with that of the other. An alternative mechanism, 
studied previously, is the valence quark-quark collision due to the hard pomeron exchange[8]. 
Both hard q-q scattering processes predict smooth decrease of da/dt in the large |t| region and 
differential cross sections which are comparable within an order of magnitude. 

We observe that the proton structure we have arrived at from our phenomenological inves- 
tigation of high energy pp and pp elastic scattering is a chiral bag embedded in a quark-antiquark 
condensed ground state. This structure (Fig. 1) is described by a nonperturbative effective field 
theory model — a gauged Gell-Mann-Levy linear sigma-model with spontaneous breakdown of 
chiral symmetry and confined valence quarks [7]. In contrast, perturbative QCD cannot provide a 
physical structure of the proton, nor can it provide quark confinement. Furthermore, the valence 
quarks of perturbative QCD are current quarks, whereas those of nonperturbative effective field 
theory models are constituent quarks. What all this implies is that perturbative QCD cannot pro- 
vide much physical insight of nonperturbative physics. As an example, the small size of the chiral 
bag R ~ 0.2 F in our case indicates that the low energy properties of the nucleon as predicted 
by this model will be essentially the same as those of the chiral topological soliton model, which 
corresponds to i? ^ 0, i.e. with no explicit quark degrees of freedom [20]. Both the chiral bag 
model and the chiral topological soliton model, for instance, will predict an axial charge distri- 
bution of size < r\ >^/^~ 0.62 F close to the experimental value [21]. From perturbative QCD, 
which has explicit quark degrees of freedom, one could not have foreseen this nonperturbative 
result. 

Finally, we note that pp elastic da/dt calculations at ^/s = 14 TeV have also been carried 
out for the whole |t| region 0-10 GeV^ using Block et al. model [17] as well as other models 
[22-24]. All these models predict visible oscillations as well as much smaller cross sections than 
ours in the large |t| region (Fig. 4). Therefore, precise measurement of elastic da/dt at large |t| 
by the TOTEM group will be able to distinguish between our model and the other models and 
shed light on the dynamics of deep-elastic pp scattering. 
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Fig. 1 . Model of proton structure underlying our study of high energy pp elastic scattering: proton has 
an outer cloud of qq condensed ground state, an inner shell of baryonic charge density and an inner core 
of valence quarks. 
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Fig. 2. Hard collision of a valence quark from one proton with a valence quark from the other proton. 




|t| (GeV^) 
Fig. 3. Solid curve shows our predicted pp elastic da/dt at LHC at y^ = 14 TeV and |i| = - 10 GeV^ 
from combined diffraction, uj exchange and a single hard valence quark-quark collision due to their low-x 
gluon clouds. Dotted curve shows da/dt due to diffraction only, dot-dashed curve that due to ui exchange 
only, and the dashed curve shows da/dt from the valence quark-quark collision only. 
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|t| (GeV^) 
Fig. 4. Comparison of our predicted da/dt (solid curve) at y/s = 14 TeV with the da/dt predicted by 
Block et al. [17], Bourrely et al. [20], Desgrolard et al. [21], and Petrov et al. (3 pomerons) [22] models. 
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